Asymmetric Zero-Bias Anomaly for Strongly Interacting Electrons in One Dimension 
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We study a system of one-dimensional electrons in the regime of strong repulsive interactions, 
where the spin exchange coupling J is small compared with the Fermi energy, and the conventional 
Tomonaga-Luttinger theory does not apply. We show that the tunneling density of states has a 
form of an asymmetric peak centered near the Fermi level. In the spin-incoherent regime, where 
the temperature is large compared to J, the density of states falls off as a power law of energy e 
measured from the Fermi level, with the prefactor at positive energies being twice as large as that 
at the negative ones. In contrast, at temperatures below J the density of states forms a split peak 
with most of the weight shifted to negative e. 
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The discovery of conductance quantization in quantum 
wires has stimulated interest in transport properties 
of one-dimensional conductors. From theoretical point of 
view, these systems are interesting because in one dimen- 
sion interacting electrons form the so-called Luttinger 
liquid 2] , with properties very different from the conven- 
tional Fermi liquids. A number of non-trivial properties 
of the Luttinger liquid, such as the power-law dependence 
of the tunneling density of states on energy and tempera- 
ture, have been recently observed experimentally 00,0. 

It is important to note that the Luttinger-liquid picture 
describes only the low-energy properties of the system. 
Quantitatively, this means that all the important energy 
scales, such as the temperature T, must be much smaller 
than the typical bandwidth of the problem. In a system 
of spin- i electrons the charge and spin excitations propa- 
gate at different velocities p , resulting in two bandwidth 
parameters. In the non-interacting case both bandwidths 



are equal to the Fermi energy Ep = (nhny 



where 



n is the electron density and m is the effective mass. Re- 
pulsive interactions between electrons increase the charge 
bandwidth D p and decrease the spin bandwidth D a . At 
moderate interaction strength both bandwidths remain 
of the order of Fermi energy Ep, and the Luttinger-liquid 
picture applies at T <C Ep. On the other hand, if the 
interactions are strong, the exchange coupling of elec- 
tron spins J is strongly suppressed, and D a ~ J -C D p . 
As a result the Luttinger-liquid picture applies only at 
T <C J, and there appears an interesting new regime 
when the temperature is in the range J <C T -C D p . In 
this regime the temperature does not strongly affect the 
charge excitations in the system, but completely destroys 
the ordering of the electron spins. 

A number of interesting phenomena have been pre- 
dicted in this so-called spin-incoherent regime. The de- 
struction of spin order may be responsible for the 
anomalous quantization of conductance in the experi- 
ments 0. Furthermore, contrary to the conventional 



Luttinger-liquid picture, the tunneling density of states 
may show a power-law peak near the Fermi level Hfl | 
even in the case of repulsive interactions, 



u{e) oc 



J < T < |e| < £> p . (1) 



Here e is the energy of the tunneling electron measured 
from the Fermi level. The result was first obtained 
for the Hubbard model in the limit of strong on- 
site repulsion U — > oo. In this case the interactions 
have a very short range, resulting in K p = 1/2 and 
v oc [|e| ln(D p /|e|)] -1 / 2 . For longer-range interactions 
K p is below 1/2, but as long as K p > 1/4, the density of 
states has a peak at low energies. 

A similar enhancement of the density of states at e <C 
Dp was predicted earlier 0] in the strongly interacting 
limit of the Hubbard model at zero temperature, 



v{e) oc 



1-3/8 



T = 0, J < lei < Dp 



(2) 



Here the effective exchange coupling of electron spins J ~ 
i 2 /U , the bandwidth D p ~ t, and t is the hopping matrix 
element in the Hubbard model. The different power-law 
behaviors of the density of states £Q) and point to the 
non-trivial physics developing when the temperature T 
is lowered below the exchange J, even if they are both 
small compared to the energy e. 

In this paper we develop a unifying theory, which en- 
ables one, in principle, to obtain the density of states in a 
system of strongly interacting one-dimensional electrons 
at arbitrary ratio T/J, as long as T, J -C e. In the spin- 
incoherent case, T> J, our theory reproduces the result 
(Q. Furthermore, we show that true asymptotic behav- 
ior of v(e) at low energies is given by Eq. (JIJ even at 
T = 0. Most importantly, in both cases the peak in v{e) 
is very asymmetric. In particular, the 3/8-power law (J2J) 
appears for short-range interactions at moderately low 
positive energies s, but never at e < 0. 
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Our approach is based on the fact that at strong inter- 
actions, when J / Dp — > 0, the Hamiltonian of the system 
of one-dimensional electrons separates into a sum of two 
contributions, H = H p + H a , describing the charge and 
spin degrees of freedom, respectively. This result was first 
obtained by Ogata and Shiba ^3 m the U/t — » oo limit 
of the Hubbard model. In the case of quantum wires at 
low electron densities a similar separation of charge and 
spin degrees of freedom was discussed in Ref.0. This sep- 
aration occurs whenever the repulsive interactions are so 
strong that even electrons with opposite spins do not oc- 
cupy the same point in space. In this case the electrons 
are well separated from each other, and their spins form 
a spin chain. The residual coupling J of the neighbor- 
ing spins is weak, J <§C Ep, and the spin excitations are 
described by the isotropic Heisenberg model 



^JSrS, 



i+i- 



(3) 



Since the Pauli principle is effectively enforced by the 
interactions even for electrons with opposite spins, the 
charge part H p of the Hamiltonian can be written in 
terms of spinless fermions (holons). In the case of the 
Hubbard model 0] the holons are non-interacting, be- 
cause two holons never occupy the same lattice site, and 
the interaction range does not extend beyond single site. 
On the other hand, if the original interaction between 
electrons has non-zero range, the holons do interact. 
Since we are only interested in the density of states at 
energies e low compared to the bandwidth D p , the holon 
Hamiltonian can be bosonized, 



//„ = / [K{d x ef 



+ K- 1 (d x j>) 2 ] dx. 



(4) 



Here u p is the speed of charge excitations and K is related 
to the standard Luttinger-liquid parameter Kj, for the 
charge modes in an interacting electron system 2] as K — 
2K p , with the factor of 2 originating from the different 
definition of the bosonic fields 4> an d 6. 

To find the tunneling density of states, one needs an 
expression for the electron creation and annihilation op- 
erators. The electron annihilation operator ipcr(x) affects 
both the charge and spin degrees of freedom: it destroys 
a holon at point x and removes a site with spin a from 
the spin chain (|3J). Building on the ideas of Refs. HI Hoi 
ITU and we write Vv(aO as 



i><j(x) 



(2 7 ra) 1 /2 



(5) 



The first factor is the bosonized form of an operator de- 
stroying a holon on the right- or left-moving branch, de- 
pending on the sign in the exponent. [Here a = Tiu p /D p 
is the short-distance cutoff; the holon Fermi momentum 
is related to the mean electron densityn as kp = 7m.] 
The operator Z^ a introduced in Ref. Illl removes site I 



with spin a from the spin chain. It is important to note 
that Eq. (0 properly accounts for the fact 0, El that 
charge modes shift the spin chain at point x by 81 — 
(j}(x)/iT with respect to its average position I = kpx/n. 

It is convenient to express the operators Zi t<T in terms 
of their Fourier transforms z a (q), where the momentum 
q is defined on a lattice and assumes values between — it 
and 7T. Then for the operator destroying a right-moving 
holon we obtain 



* dq e ik F (i+^)x 
.^ Za{q) (27^)1/2 



,i(l+i)<£0)-i6l(x) 



(G) 



In addition to the holon destruction operator e 1 ^^ 6 " 1 the 
integrand contains a factor e^ 9 ^'^. Similar factors ap- 
pear when bosonization is applied to the problem of X- 
ray-edge singularity 0, , where they represent the ef- 
fect of the core-hole potential on the electronic wavefunc- 
tions. More specifically, a core hole with the scattering 
phase shift S adds a factor e l ( 2S / 7T )<f> to the fermion oper- 
ator. Thus according to Eq. © the electron tunneling 
process that changes the momentum of the spin chain by 
q also adds a scattering phase shift for the holons 5 = q/2. 
This observation is consistent with the fact that for 
the state of momentum Q of the spin chain, the periodic 

boundary conditions for the holons acquire a phase factor 
e iQ_ 

The tunneling density of states u(e) can be computed 
as imaginary part of the electron Green's function. We 
first consider the limit of zero temperature. At |e| ^> J 
one can neglect the time dependence of the operators 
z a {q) and obtain 



d <l cf(q) 



^ 2tt r(A(g) + l) \D 



A(g) 



(J) 



for positive and negative e. Here vq = (Trhu p ) 1 and the 
exponent X(q) is given by 



- K + — 

nJ K 



- 1. 



We have also defined 



-iql 



(8) 



(9a) 



(9b) 



with averaging performed over the ground state of the 
spin chain 

At |s| <C D p the dominant contribution to the inte- 
gral in Eq. (jTJ) comes from the vicinity of its lower limit, 
q = —it. The \e\/D p — * asymptote has the form 




In 



(10) 
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FIG. 1: Sketch of the dependences c^(q) and c„{q) at zero 
temperature (solid and dashed lines, respectively), based on 
the numerical results of Ref. ITU 



This result assumes that the functions (j5J do not vanish 
at q = ±7r. 

It is important to point out that even though we so far 
assumed T = 0, the low-energy asymptote (|1(J[) agrees 
with Eq. (JIJ, rather than Eq. (J2J. Indeed, in the case of 
the Hubbard model the holons are non-interacting, the 
parameter K = 1, and the density of states behaves as 
v (e) oc l/V^j instead of the 3/8-power law 0. 

To resolve this disagreement, we first notice that at 
K = 1 our Eqs. Q and JBJ are essentially equivalent 
to Eqs. (4), (10), and (13) of Ref. HH In the notations 
of Ref. [ll| the functions are given by c+ (q) — (JV + 
l)Ccr,jv(7r— ?) andc~(g) = (JV— 1)£> - ! jv(t'"— 9) in the limit 
when the number of sites JV in the spin chain is infinite. 
These functions have been computed numerically 
and are shown schematically in Fig. ^ It is worth noting 
that the function c£(q) is numerically small at tt/2 < 
lq| < 7T, with c+(±7r) » 0.045, whereas c~(±7r) ps 0.46 
jlfij . This indicates that at very low energies e (but still 
£ ^> J) the peak of the density of states is very 
asymmetric, with the tail below the Fermi level being 
higher than the one above it by an order of magnitude. 

Given the numerical smallness of the leading asymp- 
tote (|10|l of the density of states at positive e, it is worth 
considering the subleading contributions to the integral 
in Eq. Q. They come from the values of q in the range 
— 7r/2 < q < tt/2, where c+(q) is of order one. Taking 
into account the fact that c£(q) di verg es as xl \fk~+ tt/2 
with x ~ 0.8 at q -> -tt/2 + doEi, we find 




ln^ 



(11) 

becomes —3/8, in 



At K = 1 the exponent in Eq. 
agreement with Eq. J5J. 

At low positive energies the density of states can be 



by Eqs. (jTUj) and fTTft . 



- a i given 
The subleading contribution 



diverges less rapidly than ^+ at e/D p — > 0, but with the 
numerical coefficient that is larger by a factor of about 
20. Thus for practical purposes the peak in the density 
of states v a {e) is given by v+(e), Eq. JUJ, at positive e, 
and by v^(e), Eq. (|T0|> . at negative e. 

We now turn to the spin-incoherent regime, T ^> J. 
Assuming that \e\ ^> T, one can still use Eq. (0, how- 
ever, the definitions © of the functions c^{q) should 
now assume ensemble averaging. At T ^> J the func- 
tions c^(g) can be ea sily computed by using the follow- 
ing simple argument [19j . Given that the operators Z l 

and Z\ G remove and add a site with spin a at posi- 
tion /, it is clear that the ensemble average (Zq a Z l ) 
equals the probability of all the spins on sites 0,1, ... ,1 
being a. At J -C T the spins are completely random, 
so = l/2l'l+ 1 . Similarly, {Z la Z\. a ) = 1/2 |(| . 

Then the definitions jnj) give 



ct(q) = 2c-(q) 



5 — 4 cos g 



(12) 



The expression for c„ (q) is equivalent to the result for 



3 CT falis 
.ef-Ha 



D a ,N{Q) found in Ref. 

It is important to point out that c+(q) and c~(q) differ 
by a factor of 2. As a result, the density of states l|10|) has 
a clear asymmetry: v a (e) — 2v a {—e) at T <C e <C J3 p . 
The physical meaning of this result is very simple: the 
probabilities of adding an electron with spin a at energy 
e and removing one at — e differ by a factor of 2 because 
the electron that is being removed has the correct spin 
with probability 1/2. 

The tunneling density of states can be studied experi- 
mentally by measuring the I-V characteristics of tunnel- 
ing junctions in which one of the leads is a quantum wire. 
When the electron density in the wire is sufficiently low, 
the exchange coupling is expected to be exponentially 
weak, and the regime J ~ T D p can be achieved 
0- In the experiment [2(| such a measurement was per- 
formed in a situation where the second lead is another 
quantum wire. By applying magnetic field the authors 
have been able to observe momentum-resolved tunneling. 
To measure the density of states, it is more convenient 
to make a point junction from a metal tip to the side of 
the quantum wire. Contrary to the expectations based 
on the Luttinger-liquid theory, we predict that the I-V 
characteristic of such a junction will be very asymmetric 
with respect to reversal of the bias when J <C eV <g; D p . 
The dependence v a {eV) oc dl/dV should have a peak 
at low bias at K p > 1/4, or a dip at K p < 1/4, but 
the asymmetry should be observed in either case. This 
asymmetry will indicate that electron-electron scattering 
processes have become strong enough to suppress the ex- 
change coupling J of electron spins. 

The predicted asymmetry of the density of states is 
caused by the non-trivial interplay of the spin and charge 
degrees of freedom. As a result, the asymmetry should 
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FIG. 2: Sketch of the tunneling density of states v(s) in 
two regimes: high temperature T > J (solid line), and low 
temperature T < J (dashed line). At T > J and e > T 
we predict v(e) = 2v{—e). As the temperature is lowered 
below the exchange constant J, the density of states at e < 
grows by about a factor of 3. Conversely, at e > the density 
of states decreases dramatically. At |e| <C J the standard 
Luttinger liquid effects give rise to power-law suppression of 
the density of states 0], leading to a dip at low e. 



density of states at positive energies decreases. Due to 
the importance of the subleading contribution (fTT|) the 
experimental results cannot be analyzed in terms of the 
behavior of c+(7r), Eq. (fTT)f) . However, in the limit T — > 0, 
the density of states at sufficiently low positive energies 
will become much smaller than at the negative ones. 

To summarize, we predict the density of states v(e) 
in a strongly interacting system of one-dimensional elec- 
trons to show asymmetric behavior illustrated in Fig. 
The asymmetry is strongly temperature-dependent even 
at T -C |e|. Our results can be tested in experiments 
measuring tunneling from a metal tip to the side of a 
low-density quantum wire. 
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disappear in a polarizing magnetic field, fisB 3> T, J. 
Indeed, if all spins are fixed in the j direction, one easily 
finds C| (q) = 2ir5(q). Then according to Eq. J7J) the den- 
sity of states v-\{e) ~ u Q {\e\/ D p ) x ^ . This result recovers 
the standard Luttinger liquid suppression of the density 
of states and shows no asymmetry around the Fermi 
level. 

Conductance of quantum wires is expected to depend 
strongly on the ratio J/T 0. Our theory provides a 
new way to probe this ratio by observing the asymme- 
try of the density of states v{e). The signature of the 
spin-incoherent regime, J/T <C 1, is the doubling of the 
density of states at positive energies, compared to the 
negative ones, v(e) = 2v(—e). At J/T 1 the asym- 
metry is inverted, v(—e) > v(e). This evolution of the 
density of states is described by Eq. Q, where the tem- 
perature dependence is contained in the functions (q) . 
Using Eq. (0 one can easily check that in the absence of 
magnetic field the integral of c£(q) over all q is always 
larger than the integral of c~(q) by a factor of two. At 
high temperatures c£(q) = 2c~(q), but as T becomes 
lower than J, the weight is redistributed so that c+ (q) is 
large at \q\ < ir/2, whereas c~(q) is large at |g| > 7r/2, see 
Fig. ^ Since the expression (JJJ emphasizes larger values 
of \q\, the density of states v(—e) below the Fermi level 
increases, whereas v{e) decreases. 

As the temperature is reduced, c~(n) grows from 1/6 
at T > J, Eq. to about 0.46 at T/J -4 0. 

Thus the density of states at negative energies will in- 
crease by nearly a factor of 3. At the same time the 
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